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NOTATION 


A  Attenuation 

«i  ^PVPsPl 

B  b(v)2/ 2  is  Rayloigh’s  dissipation  function 

b  Dissipation  coefficient 

C  Generalized  driving  function 

c  Velocity  of  sound 

D  Thermal  diffusivity 

d  Thickness  of  bubble  screen 

F  ■  Incident  sound  energy 

Bf  Reflected  sound  energy 

F  Characteristic  frequency  of  the  gas  bubble 

fM  Minnaert’s  resonant  frequency 

f0  Resonant  frequency 

G  Universal  gas  constant 

g  A  factor  which  takes  into  account  the  effect  of  surface  tension 

h  Height  above  the  bubble  producers 

/  Idem  factor 

K  Thermal  conductivity 

k  Restoring  stiffness 

L  Lagrangian  function 

m'  Mass  of  the  gas  contained  in  the  volume  v' 

nij  Mass  of  the  gas  in  the  bubble 

m2  Genoralizod  mnss 

n0  Average  number  of  bubbles  per  unit  volume 

p '  Complox  amplitude  of  sinusoidal  pressure  »’ 

P"  Instantaneous  pressure  on  tho  bubble  surface 

/’;  Pressure  inside  the  bubble 

P0  Static  pressure 


I\  Instantaneous  pressure  in  the  undisturbed  liquid 

p  Sinusoidal  pressure  on  the  liquid  surface 

p"  Sinusoidal  pressure  on  the  bubble  surface 

pa  Acoustic  pressure  on  the  bubble  surface 

pi  Incident  sound  pressure 

pr  Reflected  sound  pressure 

Q  Number  of  cycles  required  for  the  amplitude  of  motion  to  reduce  to  e~n  of  its 

original  value 

q  Amount  of  heat  energy  transferred 

R  Radial  distance 

R '  Nonresonant  bubble  radius 

R0  Mean  bubble  radius 

ft,  Instantaneous  bubble  radius 

r  Change  in  radius  from  the  mean  bubble  radius 

S  Net  stress  dyadic 

sp  Specific  heat  at  constant  pressure 

sv  Specific  heat  at  constant  volume 

T  Absolute  temperature 

T0  Equilibrium  absolute  temperature 

t  Time 

U  Internal  energy 

V0  Equilibrium  bubble  volume 

Vj  Instantaneous  bubble  volume 

v  Change  in  volume  from  the  equilibrium  bubble  volume 

v'  Infinitesimal  element  of  volume  in  the  gas  bubble 

11'  Work  done  on  the  bubble 

.Y  Rate  of  pure  strain  dyadic 

y  Change  in  temperature  of  the  gas  as  a  function  of  tho  bubble  radius 

7. Tube  radius 

a  A  factor  which  describes  the  departure  of  tho  bubble  stiffness  from  the  adiabatic 

stiffness 
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)-  Ratio  of  specific  heals 

5  1  ,'Q  ,  ihe  damping  constant 

60  Resonant  damping  constant 

t  Angle  between  the  incident  sound  ray  and  the  normal  to  the  bubble  screen 

t)  Polytropic  exponent 

0  Change  in  temperature  from  the  equilibrium  temperature 

A  Natural  logarithmic  decrement 

A  Wavelength 

n  Coefficient  of  viscosity 

p  Density 

a  Surface  tension 

4>  (co/2D)'i 
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V.  Velocity  potential 
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ABSTRACT 


\  theoretical  discussion  of  thermal,  radiation,  and  viscous  damping  for 
resonant  air  bubbles  in  water  is  presented.  \n  error  in  the  derivation  by  Pfriem 
for  the  thermal  damping  constant  is  corrected.  The  experimental  results  verify 
that  the  damping  constant  at  resonance  is  the  sum  of  the  thermal  and  radiation 
damping,  and  possibly  viscous  damping. 


INTRODUCTION 


The  earliest  reference  to  bubbles  as  sound  sources  was  made  by  Bragg,1  who  attributed 
to  entrained  air  bubbles  the  murmuring  of  a  brook  and  the  “plunk”  of  droplets  falling  into  water. 
Minnacrt3  has  since  shown  that  the  sound  generated  by  gas  bubbles  in  liquids  is  associated 
with  simple  volume  pulsations  of  the  bubble  without  change  of  shape.  The  bubble  behaves  as 
a  simple  damped  oscillating  system  with  one  degree  of  freedom.  Therefore,  the  differential 
equation  of  motion  for  the  hu'oblo  system  has  the  same  form  as  the  second-order  linear  differ¬ 
ential  equation  for  a  mass  fastened  to  a  spring.  As  the  bubble  periodically  expands  and  con¬ 
tracts,  the  surrounding  liquid  is  the  inert  mass  which  is  set  into  vibration,  while  the  stiffnoss 
is  due  to  the  gas  in  the  bubble.  This  zero-oider  radiator  has  a  sharply  defined  resonance  at 
the  frequency: 


i  Jzy‘'o 


[11 


where  A’0  is  the  mean  radius  of  the  bubble, 

/'0  is  the  static  pressure  at  which  the  bubble  has  the  mean  radius  /f0 , 
y  is  the  ratio  of  the  specific  heats  of  the  gas  enclosed  in  the  bubble,  and 

Pj  is  the  density  of  the  liquid.  (In  this  report  a  subscript  1  will  refor  to  properties  of 
the  gas  while  subscript  2  will  refor  to  those  properties  of  the  liquid.) 

liquation  [1]  will  be  derived  later  on  page  This  resonant  frequency  derived  by  Minnaert 
assumes  an  adiabatic  equation  of  state  for  the  gas  in  the  bubble. 

The  volume  pulsation  frequency  of  nonsphorical  gas  bubbles  in  liquids  has  been  con- 
-i  lerel  by  Strasberg,3  who  used  oblate  spheroids  to  approximate  the  nonsphorical  shapes. 
This  determination  indicates  that  the  frequency  is  only  slightly  dependent  upon  the  ratio  of 
the  major  to  the  minor  axis  of  the  spheroid.  In  fact,  for  a  ratio  of  two,  the  volume  pulsation 
frequency  of  an  oblate  spheroid  differs  by  only  2  percent  from  that  of  a  sphere  with  the  same- 
volume.  Observations  have  shown  that  large  bubbles  are  generally  nonsphorical  whereas  very 
small  hubbies  tend  to  be  spherical. 

bv,-fcr,-m*,-s  arc  listed  on  paRO  34. 


In  addition  to  simple  volume  pulsations,  there  also  may  be  oscillations  in  the  shape  of 
the  bubble.  The  natural  frequency  for  the  highermodos  of  shop*  oscillation  has  been  calculated 
by  Lamb;4  ritrasberg5  has  used  this  analysis  to  demonstrate  that  shape  oscillations  do  not  seem 
to  result  in  significant  sound  pressures  except  perhaps  very  closo  to  the  bubble.  Physically, 
the  case  of  a  quadrupole  demonstrates  the  reason  for  the  feeble  sound.  The  quadrupole  repre¬ 
sents  two  sets  of  sources  and  sinks  for  sound  on  the  bubble  surface.  The  distanco  between 
each  source  and  sink  is  very  small  compared  to  the  wavelength  of  the  sound;  therefore,  on  the 
bubble  surface,  almost  all  the  sound  from  the  source  is  fed  back  into  the  sink.  The  result  is 
that,  away  from  the  bubble,  only  insignificant  sound  pressures  occur. 

The  sound  pressure  resulting  from  excitation  of  volume  pulsations  by  several  mechanisms 
has  recently  been  discussed  in  the  literature.6  The  mechanisms,  which  cause  bubbles  to  pul¬ 
sate  and  radiate  sound,  are  bubble  formation,  coalescence,  or  division;  the  motion  of  a  free 
stream  of  liquid  containing  entrained  gas  bubbles  past  an  obstacle,  or  the  flow  of  liquid  con¬ 
taining  entrained  buhblos  through  a  pipe  past  a  constriction;  and  an  incident  sound  wave. 

Experiments  conducted  by  Sorenson7  showed  that  liquids  containing  a  gas  possess 
higher  sound  damping  characteristics  than  do  these  which  are  gas-free.  Just  a  few  widely 
dispersed  bubbles  which  are  so  small  as  to  be  invisible  can  have  an  appreciable  acoustic 
effect.  Whon  a  large  numbor  of  these  small  bubbles  are  present,  the  liquid  will  bo  nearly  opaque 
acoustically.  Small  impurities  in  liquids,  such  as  suspended  particles,  have  negligible  in¬ 
fluence  in  comparison  with  the  damping  increase  due  to  bubble^.  Therefore,  bubbles  have  a 
considerable  importance  in  the  transmission  of  underwater  sound.  In  order  to  understand  the 
attenuation  of  sound  by  gas  bubbles  in  liquids,  fundamental  processes  by  which  pulsating 
bubbles  dissipate  their  energy  must  be  known.  This  report  will  investigate  the  portion  of  the 
onorgy  radiated  in  the  form  of  spherical  sound  wavos,  the  part  which  is  transformed  into  heat 
during  the  polytropic  compressions  and  expansions  of  the  enclosed  gas,  and  the  part  of  tho 
energy  lost  in  viscous  dissipation.  It  may  be  that  these  three  processes  completely  account 
for  tho  total  damping  of  gas  bubbles  in  liquids. 

THEORY 

Poriodic  onforced  changes  in  the  pressure  on  a  bubble  result  in  volumo  pulsations  of  the 
bubble.  If  tho  amplitude  of  tho  volume  oscillation  is  small,  the  motion  of  tho  bubblo  system  is 
described  by  a  socond-order  linear  differential  equation.  For  this  svstom  possessing  ono  dogreo 
of  froodom,  the  condition  of  the  bubble  system  is  defined  by  the  change  in  volume  v  from  tho 
equilibrium  volume  F0  .  Tho  instantaneous  volume  Vl  of  the  bubble  is  tho  algebraic  sum  of  tho 
mean  volumo  F0  and  v.  In  a  similar  manner,  the  instantanoous  radius  /ij  of  tho  bubble  is  tho 
nlgobrnic  sum  of  the  mean  radius  H0  and  the  radial  increment  r.  Tho  bubblo  is  assuinod  to  be 
in  an  incompressible  liquid;  at  tho  surface  of  tho  liquid,  a  sinusoidal  prossuro  p  is  applied: 

p  =  P  eiut  [21 


o 


whore  /’  is  a  constant.  Liquids  are  slightly  compressible,  but,  as  long  ns  the  bubble  size  is 
small  compared  to  the  wavelength  of  the  pressure  wave,  the  liquid  is  considered  incompress¬ 
ible.8  The  instantaneous  pressure  /’2  in  the  undisturbed  liquid  is  the  sum  of  the  sinusoidal 
pressure  Pe,cot  and  the  static  pressure  /’0: 

P2  =  Peiat  +  P0  [3] 

Towever,  at  the  bubble  surface,  the  instantaneous  pressure  P2  is  the  instantaneous  pressure 
f  2  in  the  undisturbed  liquid  minus  the  inertial  reaction  of  the  liquid  in  motion  about  the  bubble. 
For  the  moment,  until  the  inertial  reaction  of  the  liquid  is  determined,  the  instantaneous  pros- 
-ure  /' j  at  the  bubble  surface  is  defined  as  the  sum  of  the  sinusoidal  pressure  p'and  the 
static  pressure  PQ  : 

P't-P'+Po-P'e^'  +  Po  f« 

where  /’'is  the  complex  amplitude  of  the  driving  pressure  p'.  The  bubble,  which  is  in  this 
uniform  but  alternating  pressure  field,  cannot  be  in  equilibrium  with  this  oscillating  pressure 
unless  the  bubble  itself  is  pulsating.  Uniform  pressure  in  the  gas  bubble  implies  that  the 
inertia  of  the  gas  is  negligible.  The  liquid  surrounding  the  bubble  provides  the  inertia  for  the 
bubble  system.  The  equation  of  motion  for  the  bubble  system  is  derived  in  terms  of  generalized 
coordinates  by  using  Lagrange’s  equations.  When  there  are  no  dissipation  or  forcing  pressures 
present,  Lagrange’s  equations  are  written  in  terms  of  the  Lagrangian  function  L,  which  is  de¬ 
fined  as  the  kinetic  energy  minus  the  potential  energy  of  the  system: 


d 

It 


dL 

dv 


0 


15) 


When  dissipation  is  present,  the  dissipation  pressure  is  assumed  to  be  proportional  to  the 
bubble  volume  velocity  v.  Dissipation  of  this  type  may  bo  derived  in  terms  of  a  function  B, 
known  as  itayloigh’s  dissipation  function,  and  defined  as9 


B  = 


b(v)2 


(61 


whore  b  is  the  dissipation  coefficient.  The  equation  of  motion  for  the  bubble  system  when 
there  are  dissipation  and  a  generalized  driving  function  C,  whore  neither  arises  from  a  poten¬ 
tial,  is 


♦— -C  [71 

dt  \  dv  I  dv  dv 

The  potential  energy  of  the  bubble  system  is  obtained  by  assuming  that  the  gas  in  the  bubblo 
undergoes  an  adiabatic  process  during  the  volume  pulsations  of  the  bubble: 
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/»'  t  '  -  />  r  > 
1  2  1 1  ' 0  f0 


<f  r;  - - di'j ,  or 


2  O'  ! 

0 


Theroforc,  the  potential  energy  is 


rv ,  )  /'0 

P.E.  =  -j  (t>2'.  P0)dv=  — 


\s  the  bubble  periodically  expands  and  contracts,  the  surrounding  liquid  is  set  into 
vibration.  The  maximum  kinetic  energy  of  the  liquid  particles  occurs  at  the  moment  the  bubble 
has  again  recovered  its  equilibrium  volume  1’0.  The  flow  of  the  liquid  is  irrotational:  therefore, 
a  velocity  potential  exists.  The  velocity  potential  of  a  liquid  particle  at  a  distance  A’,  due  to 
a  simple  source,  in  a  liquid  at  rest  at  infinity,  is10 


and  the  velocity  of  this  liquid  particle  is 


U  =  -  V  n  =  ■ 


The  kinetic  energy  of  all  the  liquid  volume  elements  of  density  p2  is 


''2  F  . 

K.K.  =—  I  (It)2  in  H2 


The  integration  is  extended  to  infinity  because  the  bubble  i--  assumed  to  be  surrounded  by  a 
very  large  liquid  volume.  Cpon  integration,  the  above  expression  yields  the  kinetic  energy  as 

l>2 

K-K.  =  - — —  O')2  Ilf.: 

bn  /f0 

\ceordinglv,  the  I.agrangian  /.  is 


" i  .  ,  '•  o  , 

l  --  (V)2  -  f2, 

Sir/.*  2\\ 
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and  the  equation  of  motion  for  the  bubble  .-istom,  a  hen  a  sinusoidal  pressure  /V,<J<  is  applied 
at  the  surface  cf  the  liquid  and  dissipation  is  present,  is 


■Ul'„ 


l--V„ 


«■  *  ’iv  - 


=  -  /V°' 
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The  forcing  function  P  e,,‘>t  is  preceded  by  a  minus  sign  as  a  decrease  in  pressure  results  in 
an  increase  in  the  buoble  volume.  The  term  p,  I  -  H0  is  the  generalized  mass  m-,  of  the  bubble 
system.  The  stiffness  of  the  bubble  sy stem  is  defined  as  the  change  in  pressure  on  the  bubble 
surface  associated  with  the  change  in  bubble  volume:  thereiore,  the  term  >/’0  '  1'0  is  the  adia¬ 
batic  stiffness  ,ta<): 


k 


ad 


<»'2  y'o 
dr,  »o 


(181 


Therefore,  the  linear  second-order  differential  equat.on  of  motion  for  the  bubble  system  is 
written  as 


Wjt*  +  iv  +  £ad  v  -  -  PeIOJI 


1 18l 


Ven  the  bubble  i-  slightly  noiispherical,  each  term  in  Equation  ll!)l  is  nearly  independent  of 
-'nape  when  the  mean  radius  HQ  is  taken  as  the  radius  of  a  sphere  of  the  same  volume.1 1 
Transient  volume  pulsations  are  given  by  the  solution  of  liquation  [  19 J  when  the  right  side  of 
the  equation  is  -et  equal  to  aero.  Furthermore,  if  the  dissipation  is  negligible,  Equation  1 1 91 
becomes 


*  ad 


v  =  0 


120) 


and  the  resonant  frequency  of  the  bubble  system 


is 


_L 

2  n  11^1  p2 


(211 


which  is  Minnaerl’s  exnression  as  given  in  liquation  111  on  page  1. 

When  the  right  side  of  liquation  (191  is  zero,  i.c.,  the  driving  pressure  has  been  removed, 
the  sound  pulse  from  the  bubble  consists  of  a  damped  exponential  sinusoidal  oscillation.  The 
number  of  cycles  required  for  the  amplitude  of  motion  to  reduce  to  e~n  of  its  original  value  is 
the  I,*  of  the  bubble  system.  When  the  dissipation  is  small,  the  difference  between  the  frequency 
of  the  oscillation  and  the  resonant  frequency  of  ihe  bubble  sy.-tem  without  dissipation  i.,  neg¬ 
ligible;  the  0  of  the  bubble  system  is  expressed  as 


viiere  /'0  is  tho  resonant  frequency.  The  Q  may  also  be  defined  for  a  driven  system  as 

fo 


O  r.  . 


h  ~  f\ 


(231 


where  /2  and  f  are  the  two  frequencies  respectively  above  and  below  resonance  at  which  the 
average  sound  power  of  the  bubble  has  dropped  to  one-half  its  resonance  value.12  The  total 
damping  constant  S  is  now  defined  as  the  reciprocal  Q  of  the  bubble  system  or  n~  1  times  the 
natural  logarithmic  decrement  A: 
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-I  A 

4>  =  * 


(24) 


In  this  report,  the  damping  constant  in  all  cases  will  refer  to  the  reciprocal  Q  of  the  bubble. 
Tho  total  damping  may  be  explained  by  losses  originating  from  three  processes: 

1.  Thermal  damping  Sth  due  to  the  thermal  conduction  between  the  gas  in  the  bubble  and 
the  surrounding  liquid. 

2.  Sound  radiation  damping  5rtd. 

3.  Viscous  damping  5yis  due  to  viscous  forces  at  the  gas-liquid  interface. 

The  total  damping,  expressed  as  the  sum  of  these  three  processes,  is 

^  5  =  S,h  +  Srad  ♦  Svis  (25) 


THERMAL  DAMPING 

In  the  derivation  of  Equation  [li  for  the  resonant  frequency,  the  adiabatic  equation  of 
state  was  assumed.  The  pressure  and  volume  changes  are  in  phase  with  or.o  another  so  that 
(//'2//’0  equals  -ydVx/VQ.  For  tho  adiabatic  case,  there  is  no  transfor  of  heat.  In  tho  othor 
limiting  c«se  of  a  purely  isothermal  process  in  tho  gas  space,  tho  prossure  and  volumo  changos 
aro  again  in  phase;  dP'2/Pg  oquals  -dVt/Vg.  For  this  case,  thore  is  just  as  much  heat  flow¬ 
ing  outward  from  the  bubble  during  compression  as  flows  inward  during  expansion.  Tho  work 
done  by  the  driving  pressure  in  compressing  tho  gas  space  is  just  oqual  to  tho  work  done  by 
tho  oxpanding  gas  in  moving  tho  surrounding  liquid.  However,  for  the  case  of  a  real  bubble, 
the  gas  in  contact  with  the  liquid  closely  follows  the  isothermal  equation  of  state  since  the 
liquid  has  a  largo  specific  heat  and  thermal  conductivity.  In  the  renter  of  a  real  bubble  away 
from  a  substance  having  a  high  spocific  heat,  tho  gas  ne>jdy  follows  an  adiabatic  equntion  of 
state.  Therefore,  the  thermal  process  is  polytropic  for  a  real  bubble,  and  a  phase  difference 
exists  between  the  increase  in  prossure  por  unit  original  pressure  and  the  decrease  in  volume 
per  unit  original  volumo.  This  phase  difforonco  causes  a  hysteresis  effect.  The  work  done 
on  the  gas  volumo  by  the  driving  pressure  during  compression  is  more  than  the  work  done  by 
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t ;■»'  -pace  in  i.iouni;  tho  -.urrounding  liquid  during  expansion.  Tins  difference  in  the  work 
done  represents  a  net  flow  of  heat  into  the  liquid.  The  net  flow  o'  neat  into  the  liquid  is 
characterized  by  the  thermal  damping  constant. 

The  subject  of  thermal  damping  has  been  investigated  independently  by  Pfriem,'3 
Willis,14  and  Saneyosi,1^  and  all  hate  obtained  similar  results.  The  results  of  both  Willis  and 
Saneyosi  are  available,  but  unfortunately  their  derivations  are  not  easily  accessible.  There¬ 
fore,  the  derivation  as  outlined  by  Pfriem  will  be  followed. 

In  deriving  the  expression  for  the  thermal  damping  constant,  the  gas  bubble  is  assumed 
to  be  in  an  incompressible  liquid,  and  is  excited  to  volume  pulsations  by  a  sinusoidal  pressure 
P'e,0‘‘  applied  at  the  surface  of  the  bubble.  The  liquid  has  a  large  specific  heat  and  thermal 
conductivity,  and  behaves  as  a  heat  reservoir.  This  very  large  mass  of  liquid  is  capable  of 
absorbing  or  rejecting  an  unlimited  supply  of  heat  without  suffering  appreciable  changes  in  its 
temperature.  Consequently,  in  the  liquid  adjacent  to  the  gas-liquid  interface,  it  will  bo  assumed 
that  there  are  no  changes  in  temperature.  A  temperature  field  in  the  gas  will  be  found  that 
-atisfies  this  condition,  and,  the  second  condition,  that  the  temperature  in  the  center  of  the 
bubble  is  finite.  The  oscillations  in  the  pressure,  volume,  and  temperature  of  the  gas  in  the 
bubble  will  be  assumed  small.  Consequently,  the  equations  relating  those  three  thermodynamic 
coordinates  are  linear.  In  addition,  the  density  and  the  specific  heats  of  the  gas  are  regarded 
as  constant.  In  the  gas,  the  pressure  is  not  a  function  of  position  but  only  of  time.  Therefore, 
the  gas  is  in  a  uniform  but  alternating  pressure  field;  the  inertia  of  the  gas  in  the  bubble  is 
negligible.  The  heat  transfer  process  is  conduction.  Convection  is  unimportant  as  the  time 
factor  for  establishment  of  this  process  is  considerably  larger  than  the  time  consumed  during 
a  half-cycle  compression  of  the  bubble. 

In  the  subsequent  discussion,  the  thermal  damping  constant  is  found  by  first  finding  an 
expression  for  the  change  in  bubble  volume  v.  During  a  compression  of  the  bubble  by  the 
for.  ing  pressure  P 'e,(ot  applied  at  the  surface  of  the  bubble,  work  is  done  on  an  infinitesimal 
element  of  volume  in  the  gas.  The  internal  energy  of  the  gas  in  this  volume  element  is  in¬ 
creased,  and  there  is  a  flow  of  heat  from  this  volume.  A  differential  oquation  for  conservation 
of  energy  is  formed,  and  the  solution  of  this  equation,  subject  to  boundary  conditions,  yields 
the  temperature  field  for  all  points  in  the  gas  space.  Since  the  tomperature  field  is  known, 
the  change  in  bubble  volume  v  is  calculated  for  a  change  in  the  temperature  and  a  change  in 
the  driving  pressure  at  tho  bubble  surface.  The  expression  for  the  change  in  bubble  volume 
is  then  substituted  into  the  differential  equation  of  motion  for  tho  total  gas  spaco  within  the 
bubble.  This  operation  yields  the  thermal  damping  constant  and  the  stiffnoss  of  tho  gas. 
Finally,  the  thermal  damping  constant  is  determined  as  a  function  of  the  resonant  frequenev 
of  the  bubble  system. 

In  order  to  calculate  the  thermal  contributions  to  the  total  damping  of  tho  bubble  system, 
the  change  in  bubble  volume  must  first  be  determined.  When  the  driving  pressuro  at  the  bubble 
surface  compresses  the  bubble,  work  is  done  on  the  gas  space.  This  work  done  on  the  gas 


.'■pace  increases  the  internal  energy  of  the  "as,  and  also  results  in  a  transfer  of  heat  energy 
through  the  gas.  The  added  heat  is  transferred  by  conduction  fr-  n>  the  gas  bubble  into  the 
surrounding  liquid.  The  compression  process  must  obey  the  conservation  of  energy  principle 
as  stated  in  the  first  law  of  thermodynamics: 

\  (  =  A  q  *  \  W  1261 

where  \f  is  the  increase  in  internal  energy  of  the  gas  space,  \<?  is  the  heat  added  to  the  gas 
space,  and  Ml'  is  the  work  done  on  the  gas  space.  When  each  term  in  Equation  126)  is  divided 
by  an  infinitesimal  time  \t  and  \t  is  allowed  to  approach  zero  as  a  limit,  the  rate  of  increase 
in  internal  energy  is  given  as 

dU  7Tq  7r,i’  ,  , 

_  =  (27) 

dt  dt  dt 


( \  line  drawn  through  the  differential  sign  indicates  an  inexact  differential.)  At  a  point  in  tee 
gas  space,  the  rate  at  which  work  is  done  per  unit  volume  by  the  driving  pressure  on  an  infin¬ 
itesimal  element  of  volume  v'o(  the  gas  is 


dt  v'  dt 


[28] 


'since  the  volume  decreases  during  compression,  the  term  on  tho  right  side  of  Equation  (281  is 
preceded  by  a  minus  sign  as  the  rate  of  work  done  per  unit  volume  is  positive.  For  tMs  small 
element  of  volume  at  a  jxjint,  the  rate  of  increase  in  internal  energy  per  unit  volume  is 


dU  <10 1 

rfTP>  VrfT 


[20] 


where  v|(  is  the  specific  heat  of  the  gas  at  constant  volume,  and  0l  is  the  change  in  gas  tem- 
periture  from  the  equilibrium  absolute  temperature  T0.  The  rate  of  transfer  of  heat  enorgy  per 
linn  volume  for  an  infinitesimal  volume  at  a  point  as  a  result  of  conduction  is  proportional  to 
the  divergence  of  the  temperature  gradient;  the  proportionality  constant  is  the  thermal  conduc¬ 
tivity  A t  of  the  pa  ■: 

---A',  V2«,  [30] 

tit  1  1 


Since  the  temperature  is  a  function  of  time  and  radial  distance,  i.quation  I  -TO]  is  rewritten  as 

9  a  S  ‘"‘V 

---  --  -  --  I  H) 

dt  It  n>2 


where  the  spherical  coordinate  system  originate-  at  the  center  of  the  bubbh  .  When  the  ex 
pressions  for  the  rate  of  increase  in  internal  energy  per  unit  volume  a*  a  point,  the  rate  at 
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tthuh  heat  i.-  add  Oil  por  unit  volume  .it  .1  point,  and  the  :  >te  of  work  done  prr  unit  volume  at 
.1  point  are  substituted  into  Kquation  (27),  the  differential  e<|uaticn  for  the  temperature  field 
within  the  gas  space  is  obtained: 


<*,  <’(/»,)  A',  d2(R0i)  P'2  Av. 

R  dt  "  R  ’  r'  dt 


[32] 


\  substitution  can  he  made  for  the  second  term  on  the  right  side  of  the  above  equation  by  con¬ 
sidering  the  ideal  gas  equation: 


P2v'  =  m'G(TQ  +  0{)  (33) 

v  here  '  is  the  mass  of  the  gas  contained  in  the  volume  element  v',  and  G  is  the 
gas  constant.  Kquation  [33]  may  be  differentiated  with  respect  to  time  to  yield: 

dv'  m'G  MV  ^2 
2  dt  =  R  dt  ~  V  ~ 

or 

/>'  pxG  *«*,)  dP' 

v'  dt  ~  R  ~dt  df 

As  the  universal  gas  constant  G  is  the  specific  heat  at  constant  pressure,  sp  ,  minus  the 
specific  heat  at  constant  volume  and  1 

P'2  =  P0  +  p-e>“‘,  [4] 


universal 

[34] 


Kquation  (351  may  be  written  as 


''isri  <?(«(?,) 

~~R  dT 


d(R0{) 
—Jl  dt~ 


juP'ei*1 


(36) 


ConsO(|iiently,  when  this  expression  is  introduced  into  Kquation  (32) ,  the  differentia!  equation 
becomes 


d(ROt)  h\  d2(R01)  (oU 

- - - +  j - /•>V‘ui  137) 

*  l,isP1  dR2  P  1SP) 


Kquation  1371  is  a  linear  differential  equation  describing  the  temperature  field  within  the  gas 
bubble.  The  term  A’( pisp  is  a  parameter  of  the  gas  which  Kelvin  called  the  thermal  dif- 
fusivity  I)  t  of  the  gas.  Thermal  diffusivity  is  a  measure  of  the  rate  of  heat  propogation  duo 
to  the  thermal  conduction.  A  solution  of  this  differential  equation  for  the  temperature  field 
inside  the  gas  bubble  must  satisfy  certain  boundary  conditions.  \t  the  center  of  the  bubole, 


the  change  in  temperature  ",  must  '  .>  finite,  anti  the  gradient  of  the  change  in  temperature  must 
he  /.ero,  Finally,  the  change  in  temperature  must  he  zero  at  the  c.is-liquid  interface  and  tho 
gradient  of  the  change  in  temperature  must  bo  finite.  The  solution  of  liquation  [37]  may  be 
obtained  by  several  methods.  One  method  is  to  assume  that  tho  change  in  temperature  d,  is 

0I=.VP,w'  138) 

where  y  is  a  function  of  ff  only.*  Therefore,  liquation  [37]  becomes 


„  d2  cuff/" 

;o>  (ff.y)  =  D1 -  (ffy)  +  ;  - 


dff2 


piSP, 


\  p/ossible  solution  of  Equation  [39]  is 


ffy  = 


sinh  ^  l  R 
sinh  <!>  y  ff0  . 


[39] 


[40] 


Differentiation  and  substitution  of  this  expression  in  Equation  [39]  shows  that  it  is  a  solution 
if  we  identify  a,  with 

ff0ff' 


and  i h  ,  with 


Equation  [40]  is  a  solution  for  the  temperature  field  providing  the  boundary  conditions  are 
satisfied.  When  ff  equals  the  mean  bubble  radius  ffQ ,  Equation  [40]  becomos 

ffo.V=0,  y  =  0  [41] 


and  tho  gradient  of  y  is  finite:1 


•It  is  evident  that  Equation  [38]  is  not  an  exact  solution  of  the  physical  nroblcm,  although  it  satisfies  all  the 
conditions  of  the  mathematical  equations.  The  physical  reasoning  indicates  that  the  mean  temperature  inaide  the 
bubble  must  increase  toward  the  center  of  the  bubble.  However,  according  to  Equation  (.38],  at  the  bubbie  surface, 
the  temperature  gradient  is  ( dy/dR with  a  time  average  value  of  zero,  whereas  its  time  average  value  should 
be  negative  for  a  net  outward  flow  of  heat.  This  discrepancy  cornea  about  in  neglecting  certain  second-order  termc 
m  order  to  obtain  a  differential  equation  with  linear  coefficients,  and  then  in  assuming  that  0 j  and  V  are  sinusoidal. 
The  need  for  an  increase  in  the  mean  value  of  0 j  indicates  also  the  existence  of  higher  even  harmonics.  However, 
the  subsequent  treatment  of  the  solution  of  Equation! 38|  leads  to  results  that  are  quite  good.  It  is  possible  to  go 
back  and  correct  the  equations  by  *he  method  of  successive  approximations. 
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dy 

dll 


Oj  colh  <.'  J II 0 


142] 


When  II  is  very  small,  II y  is 


r  a  * i * 

U„  "»»»•  *il:o 


Consequently,  when  II  is  zero,  y  is  finite: 


y-a. 


1  1 

R0  sinh  v'/j  R0  J 


and  the  gradient  of  y  is 


^  =  0 
dll 


143) 


[44] 


[45] 


Therefore,  the  boundary  conditions  are  satisfied;  the  temperature  field  within  tho  whole  gas 
space  is  now  known.  The  change  in  bubble  volume  v  can  now  be  determined  for  a  change  in 
the  driving  pressure  at  the  bubble  surface  and  a  change  in  tho  temperature  0y  The  total 
change  in  volume  of  the  gas  space  is  the  sum  of  the  changes  in  volume  of  all  concentric  shells 
whose  radius  is  R  and  thickness  dR.  A  shell  of  thickness  dR  has  a  volume: 


v0  =  4  nil2  dll  [46] 

In  accordance  with  the  ideal  gas  law,  two  different  states  of  the  gas  are  expressed  as 

vo 

p2  [47) 

'  0 

where  T  is  the  absolute  temperature  of  the  compressed  gas.  When  Equation  [47]  is  differen¬ 
tiated,  the  rosult  is 

vo  vo 

dv  =~  dT  -  dl‘’2  [48] 

'  o  1  0 


Since 


</T  =  (?,  =  ye>mt 


138) 
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Equation  (IS)  is  rewritten  a-; 


d  e  =  4  n  c>““ 


pa  -£*l 

L  i,  '■»  J 


The  expression  for  /<" 2y  is  obtained  from  Equation  (-10): 


r ii  sinhc'j/i’ 

Pi-V.K  i'inb  0,/V0 _ 


Therefore,  the  total  change  in  the  bubble  volume  is 

5>r  ...  . 


f  »o'"  f/f2  sinh  /’Vl*2 

e  =  4  (7  «/<“*  J  <  — ; — 7-  77  -  /f - - 

0  pi  'SP,  7oL^o  sinh  (.!>  /f  J  P 


f  P*V0  3  f.  .  I 


dli  1 501 


Equation  [511  is  further  simplified  by  noting  that 


therefore, 


m\G  mi<V“-V>  P i(s„  -O 


v*v-  f  3(v-i)r  il 


»o  ^0  yP  o  1  1/;  ?  it? 


Vi  hen  the  change  in  volume  per  unit  original  volume  (V0  -  V',)  /  1'0  is  plotted  against  the  change 
in  pressure  por  unit  original  [ires sure  (P'  -  P0)/P0  on  a  pressure-volume  graph  for  the  real 
components  of  Equation  [531,  the  area  enclosed  by  the  compression  and  expansion  curves 
represents  the  net  loss  of  energy  by  heat  conduction.  Tho  work  dono  in  compressing  tho  gas 
bubble  is  more  than  the  work  done  by  the  gas  in  expanding.  The  change  in  bubble  volume  is 
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now  known.  There  remains  now  the  task  o.‘  relating  the  change  in  the  bubble  volume  and  the 
assumed  harmonic  excitation  pressure  at  the  surface  of  the  bubble  to  tho  vibrational  properties 
of  the  gas  bubble;  i.e.,  the  stiffness  and  damping  attributes. 

In  order  to  determine  the  stiffnoss  and  damping,  the  differential  equation  of  motion  for 
the  bubble  system,  which  was  given  in  Equation  [19]  on  page  5,  is  considered: 


m^v  +  v  +  kv  =  -  P  e '0}t 


[19] 


where  6(h  is  tho  thermal  dissipation  coefficient  and  the  sinusoidal  forcing  pressure  P  e,oit 
is  applied  at  the  surface  of  tho  liquid.  Since  the  inertia  of  the  gas  in  the  bubble  is  negligible, 
Equation  [19]  can  be  rewritten  for  tho  differential  equation  of  motion  for  the  total  gas  space  in 
the  bubble: 


6th  n  +  kv  -  -  (P  el01‘  +  m2ti)=  -  P'ei<at 


[5ii 


where  PVat  is  tho  sinusoidal  excitation  pressure  at  tho  surface  of  the  bubble.  When  the 
expression  for  the  change  in  bubble  volume  v: 

V0P'e'“‘  f  ,,  n  r 

— Zp - <1  +  -^— 1 ■  (<Mocoth  * ./?„)-  1  >  [62] 

' 0  l  *\Rl  L  ]J 


is  differentiated  and  substituted  into  Equation  [54],  the  following  expression  is  obtained: 


1  V0  L  t  3  (y  -  1) 

*  +  YPo\  +  *2R2 


li0  coth  </> 


[55] 


Sinco  tho  parameter  <j)  j  has  tho  symbol  /'  under  the  square  root,  which  is  undesirable,  and  there 
is  tho  need  to  soparato  Equation  [55]  into  real  and  imaginary  components,  a  substitution  for 
</»!  is  introduced: 

<!>  j2  =  (1  +  j)2  <f>?  =  2 j<f>\  =  j  — 

ui 


<!>i  -  0  +  i)  4>i  “  (l  + 

Accordingly,  the  parameter  <t>xR0  is 

•  nr  xf,plSPif 

When  the  frequency  and  the  radius  of  tho  bubble  are  kept  constant,  tho  quantity  ^11$  varies 
as  the  square  root  of  tho  density  of  tho  gas,  or  alte-nativoly  as  the  square  root  of  the  avorage 
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:r,‘"'ir<‘  inside  the  bubble  since  the  specific  heat  at  constant  pressure  and  the  thermal  con¬ 
ductivity  are  independent  of  density: 


1  o  V  P  i  WV  f  o 
(<o,li0  are  constant) 


\nother  condition  exists  for  0,/fo  when  the  excitation  frequency  is  constant  and  the  radius  of 
the  bubble  satisfies  Kquation  [ !]:• 


for  a  resonant  bubble.  Then  for  the  resonant  case,  the  parameter  0  ,/fQ  varies  as  the  average 
pressure  inside  the  bubble: 


(<u  =  u>y  =  constant) 

Hy  introducing  the  substitution  for  j  and  noting  the  identities: 

sinh  (6 ,K0  +  ;0,/fo)  =  sinh  (0,  lt0)  cos  (0,/fo)  t  j  cosh  (0,/fo)  sin  (0,tfo) 

cosh  («^  1/^0  +  j$ilt0)  =  cosh  (0,  ff0)  cos  *  j  sinh  (0,rto)  sin  (<*,/(„), 

Kquation  [55)  becomes 

sinh  (20, /f0)  -  sin  (2 0,/?o) 


*- 

k2 


-i‘obth  V0  f  3(y  -  1)  T 

+  Mth)2  yP°l  2<ilP°  L 


cosh  (20  j rt0)  -  sin  (20, /f0) 

sinh  (20  j  /f0)  +  sin  (20, /f0) 

^cosh(20,  /?0)  -  cos  (! 


W  1  \1  | 
20, «0)  "0,/foJJ  J 


[50] 


Kvcn  though  Kquation  (5G)  is  separated  into  real  and  imaginary  terms,  the  form  is  still  not 
suitable  for  determing  the  thermal  damping  constant  at  resonance.  In  Kquation  [221: 


th 


[221 


■th  0  2 


the  thermal  damping  constant  is  given  in  terms  of  the  thermal  dissipation  coefficient  &th  , 
resonant  circular  frequency  c oQ  ,  and  the  generalized  mass  m2  .  However,  the  generalized 
mass  is  simply  the  stiffness  k  divided  by  the  square  of  the  resonant  circular  frequency. 
Therefore,  the  dimensionless  thermal  damping  constant,  which  at  resonance  is  tho  reciprocal 


M 


V  or  liple<  Ike  natural  logarithmic  decrement,  is  .U  resonance,  the  maximum  value 

of  the  thermal  damping  constant  is  about  one-tenth;  therefore,  (a  6th/A)2  is  very  small  and  can 
be  neglected  with  respect  to  unity: 


Accordingly,  Equation  (56]  becomes 


£ 

I 


."Vi  f  3(y-  l)/sinh(2<6 

'~'yP^  1 1  +  2<6,tf0  \cosh  (2^)  - 


sin  (2?i,/f0)' 
cos  (2^j/f0)y 


j  sinh  (2^j/?0)  +  sin  (2<f>1R0)  j  ^ 

i-  ; 

cosh  (2<£,tf0]- cos  (20,ffo)  </.j/f0 

2(6, sinh  (2(6, R0)  -  sin  (2<6,/?0) 

y  3(y  -  1] '  cosh  (2<6,/?0>  -  cos  (2<f>1li0)j 

Therefore,  the  dimensionless  thermal  damping  constant  is 


sinh  (2(6,/f0)  +  sin  (2(6,ff0)  j 

"A«h 

cosh  (2(6, «0]  -  cos  (20, /f0)  (6,ff0 

k 

2^,/f0  sinh  (2^6, /f0)  -  sin  (2(6, ft0) 

3(y  -  1)  cosh  (26, R0)  -  cos  (2(6,/?0) 

For  large  values  of  2<5,W0,  i.e.,  equal  to  or  greater  than  5,  the  thermal  damping  constant  is 
given  to  within  1  porcont  by 


1 

1 - 

"»th_  *lg0__8(y-l) 

A  ~  2<6,ff0  2<Aj/?0 

1  + - 

3(y  -  1) 

Large  values  of  2(6, /f0  correspond  to  large  bubbles  or  the  adiabatic  case.  The  dissipation 
arising  from  heat  conduction  vanishes.  The  thormal  damping  constant  for  vory  small  values  of 
2<6,/?0,  i.e.,  equal  to  or  less  than  2,  is  given  to  within  1  porcont  by 
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[61] 


“bth_(y  -  1)  (2<ft|/?0)2 
k  y  30 


Very  small  values  of  2<£j/?0  correspond  to  small  bubbles  or  the  isothermal  case,  and  again 
dissipation  arising  from  heat  conduction  vanishes.  The  amount  of  heat  that  flows  from  the 
gas  space  into  the  liquid  during  compression  is  just  equal  to  the  amount  that  flows  from  the 
liquid  into  tho  gas  during  expansion;  the  net  flow  of  heat  is  zero.  However,  there  is  a  transi¬ 
tion  region  between  small  and  large  values  of  2<btP0  where  the  thermal  dissipation  is  a  maxi¬ 
mum.  In  this  transition  region,  the  relation  between  the  pressure  and  the  total  gas  volume  can 
bo  oxprossod  as  Vf  equals  a  constant,  where  the  exponent  ij  varies  from  unity  to  y;  the 
state  is  polytropic.  Figure  1,  a  plot  of  the  thermal  damping  constants  bth/ k  vs  2<}>lR0  , 
clearly  illustrates  this  transition  region  between  isothermal  and  adiabatic  states. 

The  stiffness  of  gas  bubbles  is  also  important  since  the  resonant  frequency  for  the 
bubble  system  is  directly  proportional  to  the  square  root  of  the  restoring  stiffness.  By  com¬ 
paring  the  real  terms  on  both  sides  of  Equation  [56],  the  stiffness  is  expressed  as 


kV0 

rpo 


i 


3(y  -  x)/sinh  (2<f>lR0)  -  sin  (2 <£,/?„) 
2d>j/?0  \cosh  (2 -  COS  (2<£,ft0) 


)} 


[62] 


For  large  values  of  2<f>  j/?0  ,  i.e.,  for  large  bubbles  and  the  stiffness  approaching  the  adiabatic 
stiffnoss,  tho  dimensionless  stiffnoss  kVQ/yP0  is  given  to  within  one  porcont  by 


Figure  1  -  Thermal  Damping  Constant  versus 
Dimensionless  Parameter  2<f>1li0 
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Pfriem,  duo  to  an  oversight,  used  only  the  first  power  of  &>6(h.'A‘  in  obtaining  Equation  [03], 
and,  therefore,  derived  the  erroneous  result: 


[64 ! 


for  the  dimensionless  stiffness.  The  stiffness  as  given  by  Equation  [63]  is  used  to  determine 
the  thermal  damping  constant  at  resonance.  Therefore,  the  result  obtained  by  Pfriem  yiolds  an 
incorrect  thermal  damping  constant  at  resonance.  For  small  values  of  2^t/?0,  i.e.,  for  small 
bubbles  and  the  stiffness  approaching  the  isothermal  stiffness,  the  dimensionless  stiffness 
is  given  to  within  one  percent  by 


’.Vhen  the  dimensionless  stiffness  £F0  /  yP0  is  plotted  as  a  function  of  as  in  Figure  2, 

the  dimensionless  stiffness  approaches  y~l  for  small  values  of  24>l II 0  and  approaches  unity 
for  largo  values  of  2<£j liQ.  The  restoring  stiffnoss  of  the  gas  and  the  thermal  damping  constant 
are  now  known.  The  only  rei  -ining  task  is  to  introduce  the  correct  expression  for  the  stiff¬ 
ness  into  the  equation  for  the  resonant  frequency,  and  thon  determine  the  thermal  damping  con¬ 
stant  Sth  at  rosonanco. 

The  correct  expression  for  the  stiffness  will  now  be  introduced  into  the  equation  for 
tho  resonant  frequency.  Minnaert  derived  the  equation  for  the  resonant  frequency  of  pulsating 
gas  bubbles  in  liquids  by  considering  an  adiabatic  oquation  of  state: 


Figuro  2  -  Dimonsionless  Stiffness  versus 
Dimonsionloss  Parameter  26j/f0 
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However,  Che  state  is  poiytropic;  the  stiffness  constant  k  for  large  bubbles,  which  is  of  pram 
tical  importance  in  discussing  underwater  sound  transmission,  is  given  by  Equation  [63). 


A- 


3(y  -  1)  /  3(y-  1)\  r 

2 Wo  \  +  26 ,  li0  )J 


1 63a) 


where 


np0  Yp  o 


a 


y  3(y  -  1)  /  3(y  —  1)  \ 

n  =  +  26,«0  l  +  26,«0  j 


[361 


The  factor  a  describes  the  departure  of  the  bubble  stiffness  from  the  adiabatic  stiffness.  Con¬ 
sequently,  Equation  [1]  becomes 


In  the  discussion  so  far,  the  instantaneous  pressure  inside  the  bubble  has  been  considered  the 
same  as  the  instantaneous  pressure  on  the  surface  of  the  bubble.  However,  when  tho  bubble 
is  small,  the  surface  tension  pressure  increases  the  pressure  inside  tho  bubble;  consequently, 
the  instantaneous  pressure  inside  the  bubble  is  greater  than  the  instantaneous  pressure  on  the 
bubble  surface.  Smith;16  Briggs,  Johnson,  and  Mason;17  Spitzer;18  and  Robinson  and 
Buchanan19  arc  among  some  of  the  investigators  who  have  discussed  the  effect  of  surface 
tension  on  the  bubble  stiffness.  The  problem  is  to  relate  tho  pressuro  on  tho  surface  of  the 
bubble,  which  is  associated  with  the  change  in  bubble  volume,  to  the  pressure  inside  the 
bubble.  The  stiffness  A:' is  defined  as 

^  *Vl 

The  pressure  Pi  inside  the  bubble  is 

P.  =  Pi  +  2s. 

1  2  «i 


[18a) 


[68] 
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''hero  /’'  is  the*  pressure  on  the  bubble  surface,  o  is  the  surface  tension,  anil  It.  is  the  instan- 
taneou'  bubble  radius.  The  poly  tropic  equation  of  state  for  the  gas  inside  the  bubble  is 


Therefore, 


and 


w  here 


'i  7$’ -tt 

.■.JIL.Jp  MM 

avt  'Vo  njKvJ  vl  3 itxvx 

nPo(  2  a  2a  \ 

m\  {  +  vT0"  Vo^oj 

l'o  yP  o<7 


k’  = -  q  = - 

r0  ‘  V 
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Therefore,  the  correct  expression  for  the  rosonant  frequency  fQ  is 


[  (591 


[70] 
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[71a] 


[71b] 


[721 


When  the  bubbles  are  very  large,  the  stiffness  is  the  adiabatic  stiffness  and  also  surface  ten 
sion  effects  are  negligible;  consequently,  the  ratio  g/a  is  unity  and  the  resonant  frequency 
is  given  exactly  by  Minnacrt’s  equation. 

Since  damping  is  of  prime  importance  at  resonance,  the  thermal  damping  constant  will 
now  bo  determined  as  a  function  of  tho  resonant  frequency.  Tho  dimensionloss  parameter 
d>j/if0  is  written  ns 
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where 


o 

F  b - 

4 npj^j 

The  quantity  F,  for  a  given  pressure,  is  a  constant  for  tho  gas.  Since  the  parameter  F  has  the 
dimensions  of  reciprocal  time,  it  is  sometimes  called  the  characteristic  frequency  of  the  gas 
bubble.  When  Equation  (73a)  is  squared  and  the  expression  for  a  substituted,  a  quadratic 
equation  for  d>j/f0  rosults  in  terms  of  the  resonant  frequency  f0,  characteristic  frequency  F, 
and  g.  The  solution  of  this  quadratic  equation  is 


=  j  (y '  # 
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9 (y-  l)2  fo 


[74] 


Fly  substituting  the  value  of  <j> ,/?0,  as  given  by  Equation  [74],  into  Equation  [60]: 


1  <M° 

k  ~  2 *,S0 
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the  thermal  damping  constant  5lh 

<Vth 
*-  — — 
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at  resonance  is  found: 
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[75] 


As  Pfriem  used  an  incorrect  expression  for  a,  he  derived  an  incorrect  thermal  damping  con¬ 
stant;  his  expression  is 
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i/3y - 1  Rg 

Vsfy-  1)  “  f0  l  (y  r  3(y  -  1)  ”  f0 


[76] 


When  tho  following  values  are  used: 

P0  m  1  x  106  dynes/  cm2, 
o  =  75  dynos/cm, 

A',  =  5.6  x  10“5  cal/cm-soc-dog  c, 
y  =  1.40 


20 


s  0.2  1  cal  gm. 


I' 2 
<>  1 


1.00  gin  cm3. 


1.2!)  •  10"  3  gn  cm3. 


the  thermal  lamping  constant  <>(h  calculated  f 
u-i:ig  liquation  [75|  and  that  using  I’friem’s  i 
equation  are  plotted  in  Figure  :i  a-  a  function  t 
of  the  resonant  frequency.  Figure  3  dearly  * 
shows  that  the  result  obtained  hv  Pfrieni  is 
50  to  (55  percent  too  high.  The  thermal  damping 
constant  determined  in  this  report,  using  the 
method  outlined  by  Pfrieni,  agrees  exactly  with 
the  curve  of  Viillis  as  given  in  the  report  by 
spitzer.  The  thermal  damping  constant  8(h  for 

air  bubbles  in  water  larger  than  15  microns  radius  and  with  resonant  frequencies  less  than  24*! 
kilocycles  |>er  second  is  determined  to  within  1  |>ercent  by  using  Equation  [75].  When  the 
resonant  frequency  is  small,  Kquation  [75]  for  the  thermal  damping  constant  is  replaced  by  the 
simple  equation: 


i 

JO 


fffQuCnCy  m  fcitOCy Clfl  per  Second 


Figure  55  -  Thermal  Damping  Constant  for 
Resonant  Air  Bubbles  in  Water 

The  upper  curve  is  the  incorrect  theoretical 
result  obtained  by  Pfricm. 


th 


g  4.41  X  10-«  V/0  (seconds)  [77] 


where  ft)  is  the  resonant  frequency  in  cycles  per  second.  This  simple  relationship  gives  the 
thermal  damping  constant  to  within  1  percent  for  nir  bubbles  in  water  larger  than  0.05  cm; 
i.o. ,  for  air  bubbles  with  resonant  frequencies  less  than  about  7  kilocycles  per  second.  Kqua- 
tion  [771  icveals  that  the  thermal  damping  constant  at  low  frequencies  is  proportional  to  the 
square  root  of  the  resonant  frequency. 

RADIATION  DAMPING 

In  a  compressible  liquid,  a  bubble  excited  intovolumo  pulsations  oxpends  a  portion  of 
its  energy  by  radiating  spherical  sound  waves.  The  bubble  is  considered  as  a  simple  sound 
source;  the  bubble  radius  ltQ  is  considered  small  compared  to  the  wavelength  \  of  tho  rudiatod 
sound.  Smith20  has  calculated  tho  radiation  damping  for  gas  bubbles  in  liquids.  In  order  to 
derive  (he  expression  for  the  radiation  damping  constant,  the  velocity  potential  for  a  simplo 
sinusoidal  source  in  a  compressible  liquid  is  stated;. 

(i)V. 

O  =; - c,<u(r  -  K/c2)  [78] 

0  Unit 
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where  c,  is  the  velocity  of  sound  in  the  liquid, 

1C  is  the  radial  distance,  and 

t’j  is  the  complex  amplitude  of  the  change  in  bubble  volume  v. 
Equation  1 52]: 


[52] 
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The  acoustic  pressure  !s  determined  from  the  velocity  potential  by  the  following  equation: 


P2(°2vi 

Fa  ~f’2  ~dl=~  4 nil 
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[81] 


On  the  bubble  surface,  the  acoustic  pressure  is 


,  -.-P^—Lei^e~,<oR°/C2 

P°  4W7ft 


[82] 


or 


Pa 


p20>2  V} 

~i TT 


[82a] 


where  only  the  first  four  tenns  are  kept  in  the  expansion.  The  acoustic  pressure  on  the  bubble 
surface  is  just  the  difference  between  the  driving  pressure  on  the  surface  of  the  linuid  and 
the  change  in  pressure  on  the  bubble  surface  associated  with  the  change  in  bubble  volume. 

The  change  in  pressure  on  the  bubble  surface  is  k'v.  Therefore, 


Pa  f  ■  *  v  -  Pa  +  k'vle,<ot  =  -  Pc>Cit 


oo 


j2  ;'j  e1  wl  j  (o2Hq\  p2<o2  vl  el<x  1  I<oR0 
4  nli0  \  c2oi/  *  4>r/f0  ^  c2 

—  ♦*>,< 

I  nil 


[841 


,e>“1  =  -  Pei 0,1 
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then 
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For  n  large  resonant  air  bubble  in  water  with  a  resonant  frequency  of  one  kilocyclo  per  second 
and  a  radius  of  0.33  cm,  the  term  a>0/l0/c2  is 

a>0lt0  2*rS0  1 

c2  \2  75 

Therefore,  the  terms  of  higher  order  than  a>0R0/  c2  will  be  neglected  and  Equation  [87]  becomes 


?2  P2 w 

v  +  - -  v  +  k'v  =  -  Pel<ot 


4a/f0  4hc2 


[88] 


\s  long  as  the  bubble  radius  is  small  compared  to  the  wavelength  of  the  radiated  sound,  it  is 
seen  that  the  generalized  mass  term  for  the  case  of  a  compressible  liquid  is  the  same  as  the 
corresponding  term  in  Equation  [17]  where  the  liquid  was  considered  incompressible-. 


P  2 

-  v  +  bv  4  k'v  =  -  Pe,cot 
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The  radiation  dissipation  coefficient  6faJ  is 


P  2 

4  rr  c2 


[17a] 


[80] 


23 


Therefore,  ihe  radiation  damping  constant  at  resonance  is 

1  Jrad  -"/l>o^o 

rad  ^rad  <e0«2  C2 

where  «2  is  the  generalized  mass.  The  factor  tj  takes  into  account  tne  effect  of  surface  ten¬ 
sion  whereas  a  describes  the  departure  of  the  bubble  stiffness  from  the  adiabatic  stiffness, 
for  large  bubbles,  the  ratio  tj.'a  is  unity,  and,  since  the  quantity  2n  ltQ  /„ c2  is  a  constant, 
the  radiation  damping  constant  at  resonance  is  independent  of  frequency. 

VISCOUS  DAMPING 

The  problem  of  a  pulsating  spherical  bubble  in  a  viscous,  incompressible  liquid  will 
now  be  treated;  the  viscous  damping  constant  at  resonance  will  be  derived.  Mallock21  in  1910, 
and,  later,  Spitzer22  and  Poritsky23  investigated  this  problem,  for  a  pulsating  bubble,  the 
effect  of  viscosity  is  perhaps  difficult  to  visualize.  Lamb24  states,  “The  only  condition  under 
which  a  liquid  can  be  in  motion  without  dissipation  of  energy  by  viscosity  is  that  there  must 
he  nowhere  any  extensions  or  contraction  of  linear  elements;  in  other  words,  the  motion  must 
consist  of  a  translation  and  a  rotation  of  the  mass  as  a  whole,  as  in  the  case  of  a  rigid  body.” 
Mallock  gives  us  a  physical  picture  of  the  effect  of  viscosity  on  a  pulsating  bubble  in  an  in¬ 
compressible,  viscous  liquid  by  considering  a  small  element  of  a  spherical  shell  of  liquid  at 
the  bubble  surface.  This  element  has  definite  radial  and  lateral  dimensions  at  the  instant  the 
bubble  radius  is  at  its  mean  position.  When  the  bubble  expands,  the  small  'liquid  clement  is 
distorted;  the  radial  thickness  decreases  while  the  lateral  dimension  increases.  Likewise, 
when  the  bubble  contracts,  the  liquid  element  is  again  distorted;  this  time  the  radial  thickness 
increases  and  the  lateral  dimension  decreases.  Since  the  liquid  is  incompressible,  the  distor¬ 
tion  is  not  caused  by  a  change  in  the  volume  of  the  liquid  element  but  by  viscous  stresses. 
Consequently,  more  energy  is  required  to  compress  the  bubble  than  is  regained  in  the  subse¬ 
quent  expansion. 

In  the  presence  of  viscosity,  momentum  is  transmitted  from  one  region  of  the  liquid  to 
another  moving  at  a  different  velocity.  An  element  of  liquid  moving  rapidly  in  a  particular 
direction  tends  to  transmit  its  momentum  to  other  elements  of  the  liquid.  The  Navier-Stokes 
equation  of  motion  describes  the  force  per  unit  volume  acting  on  an  infinitesimal  element  of 
volume,  at  a  point  in  a  viscous  liquid.  The  force  per  unit  volume  is  duo  to  the  instantaneous 
pressure  distribution  of  the  surrounding  liquid  as  in  the  case  of  a  nonviscous  liquid,  and  is 
also  due  to  the  rate  of  change  of  momentum  caused  by  the  presence  of  viscosity.  When  there 
are  no  external  forces  acting  on  each  unit  mass  of  the  liquid,  the  equation  of  motion  is 


e. 


V? 


190) 


where 

It  =  radial  velocity  vector, 

(  a  \  r»- 

—  *  K  •  V]  H  =  acceleration  vector, 

\')t  I 

p0  =  mean  pressure,  and 
H  =  coefficient  of  viscosity. 

rV 

\s  the  liquid  is  considered  to  be  incompressible,  the  divergence  of  the  velocity  V ■  H  vanishes 
so  that  the  second  tern  on  the  right  side  of  Equation  [91]  disappears.  The  only  remaining 
viscous  term  is  a  V2fc .  Since  the  motion  of  the  liquid  is  irrotational,  the  velocity  can  be  ex¬ 
pressed  as  the  gradient  of  a  scalar  velocity  potential  ft,  and  the  last  term  on  the  right  side  of 
Equation  (91)  is  written: 


p 


V 2 


it  =nv2(-v ft)  =  -p7(r.  pft) 


pv2  ft  =  nv  (v  ■  ft)  =  o 

Therefore,  there  are  no  net  viscous  forces  acting  inside  the  liquid  for  the  case  of  a  pulsating 
spherical  bubble  in  an  incompressible,  viscous  liquid.  Due  to  the  presence  of  viscosity, 
momentum  is  transmitted  through  the  liquid,  but  each  infinitesimal  element  of  liquid  volume 
receives  just  as  much  as  it  loses;  therefore,  there  is  no  net  viscous  force  acting  on  any  ele¬ 
ment  of  volume  internal  to  the  liquid.  The  Navier-Stokes  equation  is  not  applicable  for  dis¬ 
cussing  the  effect  of  viscosity  for  pulsating  spherical  bubbles. 

However,  even  though  the  net  viscous  forces  in  the  liquid  vanish,  there  are  viscous 
forces  acting  at  the  surface  of  the  bubble  whore  they  exert  an  excess  pressure.  The  net  stress 
dyadic  S  is  written  as25 

o  -> 

S  =  -~  pV-ItU  2p ,V  [!>2| 

3 

where  /  is  the  idem  factor  and  X  is  the  rate  of  pure  strain  dyadic.  (The  algebraic  signs  conform 
with  the  practice  in  elasticity  of  denoting  a  tension  as  positive  and  r>  pressure  as  negative.) 

\s  the  liquid  is  incompressible,  the  divergence  of  the  velocity  is  zero,  and  Equation  [ 92] 
simplifies  to: 

•S'  =  2  g  A  I1':)  I 

Due  to  the  radial  motion  and  spherical  symmetry,  the  principal  direction-  of  -tress  anil  rate  of 
strain  must  be  radial:  this  will  be  chosen  to  correspond  to  SR  and  .V^,  respectively: 


I  <•  r:.n*  of  pure  strain  \  g  is  the  gradient  of  the  radial  velocity. 
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lierefore,  the  radial  stress  Sfi  at  the  bubble  surface  is 


itv 
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ind  the  equation  of  motion  for  the  bubble  system,  when  the  effect  of  viscosity  is  included,  is 


m,  if  +  — —  «  +  /•')>=-  P  el  <J  * 
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and 
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Therefore,  the  viscous  damping  constant  at  resonance  is 


1  ^vis  87/10  87/1  .fa 

^vis  "qTJj  3yf'0ji  0  3y/p  W\  jt 


[08] 


whore  y  is  the  factor  which  takes  into  account  the  effect  of  surface  tension,  and  a  is  a  factor 
which  describes  the  departure  of  tho  bubble  stiffness  from  the  adiabatic  stiffneso.  The  viscous 
damping  constant  8vis  is  directly  proportional  to  tho  resonant  frequency  f0. 

The  effect  of  viscosity  is  realized  only  ti  rough  the  boundary  condition  at  tho  surface 
of  the  bubble,  rather  than  through  the  Navier-Stokcs  equation  where  tho  resultant  of  tho  vis¬ 
cosity  stresses  per  unit  volume  at  any  point  internal  to  the  liquid  vanishes. 

TOTAL  DAMPING  CONSTANT  AT  RESONANCE 

The  total  damping  for  resonant  air  bubbles  in  water  may  be  explained  by  lossos  origi¬ 
nating  from  the  energy  which  is  transformed  into  heat  during  the  polytrop.c  compressions  and 
expansions  of  the  air  in  the  bubble,  the  energy  radiated  in  the  form  of  spherical  sound  waves, 


and  the  energy  lost  in  viscous  dissipation.  The  total  damping  constant  80  at  resonance  is 
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Figure  4  is  a  plot  of  the  thermal,  radiation, 
viscous,  and  total  damping  constants  as  a 
function  of  the  resonant  frequency  fQ  for  air 
hubbies  in  water.  Damping  constants  for  air 
bubbles  in  water  ranging  in  radius  from  3 
microns  to  3  millimeters  are  displayed  in 
Figure  4.  The  thermal  damping  constant  5(h 
reaches  a  maximum  around  200  kilocycles  per 
second  while  the  radiation  damping  constant 
5rad  is  nearly  constant  up  to  600  kilocyclos 
l«r  second  and  then  very  slowly  begins  to  in¬ 
crease.  The  viscous  damping  constant  5vis  is 
of  the  same  order  of  magnitude  as  the  radiation 
damping  constant  at  180  kilocyclos  por  socond; 
just  above  1000  kilocycles  per  second,  tbo 
viscous  damping  constant  and  the  thermal 
damping  constant  have  tho  same  value. 


Resonant  Frequency  1Q  in  Kilocycles  per  second 

Figure  4  -  Theoretical  Thermal,  Radiation, 
Viscous,  and  Total  Damping  Constants 
for  Resonant  Air  Bubblos  in  Water 


EXPERIMENTAL  METHODS 

There  are  essentially  four  mothods  by  which  tho  damping  constant  §0  can  bo  determined 
experimentally.  Most  of  those  methods  are  indirect  ones  involving  the  calculation  of  the  damp¬ 
ing  constant  from  certain  measured  acoustical  properties  of  tho  bubbles. 

SUCCESSIVE  OSCILLATIONS26 

The  method  of  successive  oscillations  is  a  direct  process  for  determining  the  damping 
constant.  The  signal  from  a  hydrophone,  which  is  placed  close  to  an  oscillating  bubble,  i*> 
amplified  and  applied  to  the  input  terminals  of  acathodo-ray  oscilloscope;  tho  bubble  pulse 
appears  on  the  screon  as  a  damped  sino  wavo.  If  the  amplitude  of  successive  oscillation  is 
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plotted  on  -.emi-log  jja fx*r  as  a  function  of  the  cycle  number  of  oscillation,  the  logarithmic  de- 
crement,  ami,  therefore,  the  reciprocal  Q  can  be  determine.!  fro'  the  slope.  The  resonant  fre¬ 
quency  of  the  pulsating  bubble  is  determined  by  applying  a  signal  of  known  frequency  to  the 
input  of  the  oscilloscope  and  measuring  the  time  scale  across  the  screen. 

WIDTH  OF  THE  RESONANCE  RESPONSE27,28,29 

On  page  5,  the  Q  of  the  bubble  system  was  defined  in  terms  of  the  resonant  frequency 
and  the  two  frequencies  above  and  below  at  which  the  average  sound  power  of  the  bubble  had 
decreased  to  one-half  its  resonance  value.  Since  the  power  is  proportional  to  the  square  of 
either  the  radial  velocity  or  the  radial  displacement,  the  damping  constant  S0  can  be  found  by 
plotting  the  square  of  either  of  these  parameters  as  a  function  of  the  frequency.  The  amplitude 
of  oscillation  of  large  gas  bubbles  can  be  found  using  a  photoelectric  method,  and,  for  small 
bubbles,  the  radial  velocity  can  be  measured  using  a  kind  of  velocity-ribbon  microphone. 

1.  Photoelectric  Method.  \  single  gas  bubble  oscillating  to  a  sonic  excitation  is  illumi¬ 
nated  optically  and  the  scattered  light  measured  by  a  photoelectric  cell.  The  change  in  cross 
section  of  the  bubble  image  modulates  the  quantity  of  light  received  at  the  photocell.  The 
alternating  current  generated  in  the  circuit  of  the  photocell,  by  the  change  in  the  bubble  cross 
section,  is  amplified  and  recorded  on  a  suitable  recorder.  Ry  varying  the  sonic  excitation 
frequency  and  noting  the  changes  in  the  bubble  cross  section,  the  band  width  and  the  resonant 
frequency  can  be  determined. 

2.  'tibbon  Microphone  Method.  \  single  bubble  is  caught  on  a  small  wax  sphere  fastened 
to  a  platinum  thread  which  is  placed  between  the  poles  of  an  electromagnet.  Hie  bubble  oscil¬ 
lations  are  produced  by  a  constant  frequency  magnetostriction  projector.  \s  the  bubble  oscil¬ 
lates,  the  platinum  thread  is  carried  along  with  the  oscillations  and  this  motion  of  the  platinum 
thread  produces  an  alternating  KMF  which  is  proportional  to  the  radial  velocity  of  the  gas 
bubble.  This  arrangement  represents  a  sort  of  ribbon  microphone.  Since  the  generalized  mass 
is  greater  than  the  vibrating  part  of  the  mass  of  the  thread  or  the  wax,  the  platinum  wire  and 
the  wax  are  assumed  to  exert  negligible  influences  or  the  resonant  frequency.  The  bubble  is 
allowed  to  grow  slowly  and  its  diameter  measured  with  a  microscope;  the  voltage  produced  by 
the  ribbon  microphone  traverses  a  maximum  as  tho  diameter  of  the  bubble  increases.  There¬ 
fore,  the  resonant  frequency  and  damping  constant  can  be  determined. 

STANDING-WAVE  RATIOS  -10-31-32 

\  single  gas  bubble  is  allowed  to  rise  freely  in  u  liquid-filled  tube  and  oscillate  under 
the  influence  of  a  plane  progressive  sound  wave.  The  diameter  of  the  tube  is  less  than  half 
a  wavelength  so  that  the  sound  pressure  is  constant  over  tho  cross  section  of  tho  tube.  Dis¬ 
turbances  by  reflection  of  this  sound  wave  from  the  free  surface  are  prevented  by  an  absorp¬ 
tion  device  or  by  using  a  pulse  technique.  The  sound  energy  £(,  which  is  radiated  by  a 
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•r  tp.-iiucpr  at  the  lower  end  of  the  tube,  is  partly  reflected  In  the  1  ubble  h’r  and  recorded  In 
a  probe  h\  Irophone  arranged  between  the  transducer  and  the  bubble.  The  damping  constant 
0  can  be  measured  from  the  relative  reflection  coefficient  (Er  of  a  bubble  oscillating 

at  its  resonant  frequency: 


where  Z  is  the  radius  of  the  tube,  and  p(  and  pr  are  the  sound  pressures  of  the  incident  and 
reflected  waves,  respectively.  However,  the  energy  of  the  reflected  signal  cannot  be  referred 
to  that  of  the  direct  signal  because  of  friction  losses  occuring  in  the  tube.  In  order  to  obtain 
the  corresponding  energy  of  the  direct  signal,  the  energy  of  the  reflected  pulse  is  measured 
when  the  bubble  is  replaced  by  an  obstacle  covering  the  entire  cross  section  of  the  tube,  a 
large  bubble,  for  instance.  Therefore,  the  damping  constant  can  be  measured  from  the  standii.g- 
wave  ratio,  and  the  resonant  frequency  is  that  frequency  of  the  plane  progressive  sound  wave 
which  producos  the  maximum  oscillation  for  the  bubble. 

RESONANCE  ABSORPTION33'34 

This  method  of  determining  the  damping  constant  depends  upon  measuring  the  attenuation 
nf  sound  by  a  screen  of  bubbles.  For  bubbles  of  a  single  size,  P  »  attenuation  through  a  bubble 
•  croon  is  a  maximum  at  the  resonant  frequency  of  the  bubbles.  A  projector  and  transmission 
1  ydrophone  are  located  on  opposite  sides  of  the  bubble  screen.  Each  instrument  is  faced  to- 
a  ird  the  and  the  line  joining  them,  at  the  point  of  intersection  with  the  bubble  screen, 

f  r  -  in  angle  c  with  the  normal  to  the  screen.  In  order  to  obtain  data  as  to  the  distribution 
of  bubble-  u  cording  to  size,  the  rate  of  rise  of  bubbles,  which  is  a  function  of  the  bubble 
radiu- ,  l-  determined.  If  only  a  very  short  burst  of  bubbles  is  allowed  to  escape  from  the  bub¬ 
ble  producer-  and  the  resultant  screen  is  observed  at  a  height  h  and  time  l  Inter,  the  screen 
mtains  oniy  those  bubbles  whose  rate  of  rise  is  h/t.  When  the  bubbles  are  allowed  to  rise 
freely  but  in  these  definite  bursts  or  pulses,  attenuation  measurements  versus  time  elapsed 
ifter  (he  initiation  of  the  pulse  bubble  screen  are  made  using  the  transmission  hydrophone. 

This  method  is  repeated  for  several  projector  frequencies.  Carstensen  and  Foldy35  give  the 
resonant  damping  constant  50  in  terms  of  the  attenuation  A: 
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where  is  the  average  number  of  bubbles  per  unit  \olume, 
a  is  the  thickness  of  the  bubble  screen, 

R'  is  the  off-resonant  bubble  radius,  and 


in  deriving  liquation  [101],  Carstensen  and  Foldy  made  some  assumptions  about  the  off- 
resonance  behavior  of  the  damping  constant.  They  assumed,  for  a  bubble  screen  containing 
bubbles  of  essentially  uniform  size,  that  the  off-resonance  damping  constant  equals  (H0,'lt')380 
where  50  is  the  resonant  damping  constant.  The  distribution  in  size,  space,  and  number  of 
the  hubbies  must  be  known  to  determine  the  resonant  damping  constant. 

COMPARISON  OF  THEORY  AND  EXPERIMENTAL  RESULTS 

The  theoretical  damping  constant  and  the  experimental  values  for  the  damping  constants 
are  plotted  as  a  function  of  the  resonant  frequency  in  Figure  5.  Therefore,  Figure  5  gives  an 
indication  as  to  how  well  the  experimental  results  agree  with  the  theory  of  damping. 

Meyei  and  Tamm36  have  used  the  width  of  the  resonance  response  method  to  obtnin  the 
damping  constant;  these  results  are  extremely  high.  This  high  damping  constant  may  bo  due 
to  the  particular  conditions  of  the  experiment.  In  using  the  ribbon  microphone,  considerable 
damping  may  have  been  due  to  the  oscillation  of  the  platinum  thread  in  the  magnetic  field. 

In  addition,  the  experimenters  themselves  state  the  bubbles  appeared  dull  and  blurred  near 
the  resonance  point;  consequently,  the  diameter.-,  of  the  bubbles  could  not  have  been  measured 
accurately  with  a  microscope,  which  would  affect  the  determination  of  the  resonant  frequency. 
The  damping  constant  for  large  bubbles  was  determined  using  the  photoelectric  mothod.  For 
this  procedure,  Meyer  and  Tamm,  and  later  I.auer,37  used  a  thin  wire  annulus  to  hold  the  bub¬ 
bles  and  prevent  them  from  rising  to  the  surface  while  the  measurements  wero  being  made. 
Indeed,  the  high  damping  constants  found  by  Meyer  and  Tamm  may  be  duo  to  the  wire  annuli 
adding  to  the  damping  of  the  bubble  system.  At  low  frequencies,  the  dnmping  constants 
measured  by  I.auer  are  about  25  percent  higher  than  the  theoretical  prediction. 

Hauer,  formerly  of  the  David  Taylor  Model  Basin  and  now  at  Reeves  Instrument  Com¬ 
pany,  used  the  successive  oscillation  method  for  determining  the  dnmping  constant.  In  this 
experiment,  the  damping  constant  for  a  free  bubble  was  measured;  therefore,  there  is  no 
additional  damping  due  to  a  bubble  holder.  The  bubble  was  formed  at  a  nozzle;  the  volume 
pulsation.-  started  just  as  the  bubble  closed  and  separated  from  the  nozzle.  The  unpublished 
damping  constant  measurements  of  Bauer  are  about  20  percent  higher  than  the  theoretical 
prediction. 


Figure  5  -  Theoretical  and  Experimental  Values  of  the  Damping  Constant  for 
Resonant  Air  Rubbles  in  Water 

Points  are  from  faired  curves  through  the  experimenter's  data. 
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The  method  of  standing-wave  ratios  was  used  by  Exner,38  Exner  and  llampe,39  and 
Iiaoske1‘,°  to  determine  the  resonant  damping  constant;  the  results  of  Exner,  and  Exner  and 
llampe  agree  very  well  with  the  theoretical  curve.  According  to  the  theory,  the  viscous  damp¬ 
ing  bocomos  important  around  200  kilocycles  per  second;  Haeske  has  measured  the  damping 
constant  i"  this  frequency  range.  At  resonant  frequencies  of  200  and  300  kilocyclos  por  sec¬ 
ond,  tho  damping  constants  determined  by  Hneske  are  4  to  8  percent  lower  than  the  theoretical 
curve.  When  tho  theoretical  damping  constant  curve  does  not  include  the  viscous  damping  con¬ 
stant,  but  only  the  thermal  and  radiation  damping  constants,  the  experimental  results  of  Haeske 
are  4  to  8  percent  higher  than  the  theoretical  curve.  However,  the  measurements  by  Haeske 
are  only  accurate  to  within  10  percent.  Therefore,  a  definite  conclusluii  cannot  bo  formed  as 
to  whether  viscous  damping  contributes  or  does  not  contribute  the  total  damping.  Also,  the 
value  for  the  coefficient  of  viscosity,  which  is  used  in  determining  tho  theoretical  viscous 
damping  constant,  was  obtained  experimentally  for  steady  flow.  At  high  frequencies,  the  value 
for  tho  coefficient  of  viscosity  may  bo  considerably  snmllor  than  for  tho  steady-flow  case. 

This  subject  is  now  boing  investigated.  Some  additional  damping  experiments  in  thi«  frequency 
range  using  a  different  experimental  method  may  also  docide  this  dilemma.  Above  4  0  kilo¬ 
cycles  per  second,  Exner  an-1  llampe  very  often  found  “anomalous”  bubbles  with  much  lower 
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damping  constants  than  tho  regular  bubble*.  The  measured  resonant  frequency  did  not  agree 
with  the  frequency  calculated  from  the  measured  diameter  of  the  '  jhMe  when  Equation  |7d]. 


was  used.  The  “anomalous”  bubbles  have  higher  frequencies  than  this  equation  predicts. 

It  was  noted  that  in  almost  all  cases  the  “anomalous”  bubbles  ha  1  dust  particles  on  their 
surfaces.  This  increaso  in  resonant  frequency  could  not  be  explained  by  a  decrease  in  the 
generalized  mass  as  the  dust  particles  would  add  to  this  mass,  and  there  does  not  seem  to  be 
a  logical  explanation  for  a  possible  increase  in  the  stiffness.  Strasberg41  tentatively  sugges¬ 
ted  that  this  behavior  may  perhaps  be  associated  with  surface  oscillations  of  the  bubble  since 
for  very  small  bubbles  the  frequency  of  surface  oscillations  may  be  of  the  same  order  as  the 
frequency  of  ordinary  volume  pulsations.  The  excitation  of  surface  oscillations  by  sonic 
excitation  may  require  some  nonsyn.metry  supplied  by  the  dust  particles.  The  damping  asso¬ 
ciated  with  surface  oscillations  is  not  the  same  as  tho  damping  associated  with  volumo  pulsa¬ 
tions,  and  this  would  account  for  the  different  damping  constant  measured  for  “anomalous” 
bubbles.  When  Haeske  performed  his  experiment,  he  took  extreme  care  to  obtain  clean  exper¬ 
imental  conditions,  and  found  no  trace  of  “anomalous”  bubbles  in  the  100-300  kilocycles  per 
second  range. 

Carstensen  and  Foldy42  used  the  resonance  absorption  method  to  determine  the  damping 
constant.  The  damping  constant  results  are  very  high.  The  authors  admit  they  have  only  a 
small  amount  of  evidence  to  indicate  that  the  off-resonance  damping  constant  equals(/<’o/ft')'!<>0. 
The  off-resonance  damping  constant  was  used  in  deriving  an  expression  for  tho  resonant  damp¬ 
ing  constant  in  terms  of  the  attenuation.  In  this  mothod,  a  large  number  of  bubbles  are  presont, 
and  the  exact  distribution  in  size,  space,  and  number  is  difficult  to  determine.  Also,  there  may 
be  interaction  between  individual  bubbles;  these  interactions  arc  surely  complex  and  difficult 
to  determine.  Tho  curve  for  the  resonance  damping  constant  as  a  function  of  the  resonant  fre¬ 
quency  is  certainly  broadened  by  the  presence  of  bubbles  having  different  resonant  frequencies. 
However,  this  broadening  is  difficult  to  calculate  since  the  exact  bubble  distribution  in  size, 
space,  and  number  is  not  known  to  a  sufficient  degree. 

Excluding  the  results  obtained  by  Meyer  and  Tamm,  and  Carstensen  and  Foldy,  tho 
experimental  damping  constants  agree  very  well  with  the  theoretical  curve.  Tho  damping  con¬ 
stants  at  low  resonant  frequencies  obtained  by  flauer  and  Lauer,  using  different  experimental 
methods,  agree  quite  well  with  each  other,  but  their  results  are  higher  than  the  theory  predicts. 
\s  mentioned  on  page  1,  the  ri  sonant  frequency  is  only  slightly  different  for  nonsphorical 
bubbles  .  However,  the  viscous  damping  theory  assumed  spherical  bubbles  where  the  motion 
is  completely  radial.  Perhaps  there  is  considerably  more  viscous  damping  when  the  bubble  is 
slightly  nonspheric.nl.  In  addition  to  volume  pulsations,  there  are  also  oscillations  in  the 


shape  of  the  bubble.  itecent  work*  indicates  that  for  largo  amplitude  radial  pulsations,  there 
is  some  coupling  between  tho  radial  motion  and  the  shape  oscillation  This  may  result  in  tho 
removal  of  some  of  the  energy  associated  with  tho  pulsation.  At  high  frequencies,  the  exper¬ 
imental  data  of  Haeske  do  not  confirm  whether  or  not  viscous  damping  is  important. 

SUMMARY  AND  CONCLUSION 

Bubblos  oxcitod  to  volume  pulsations  have  a  polytropic  oquation  of  state  for  the  gas 
which  results  in  a  phase  difference  between  the  change  in  pressure  per  unit  original  pressure 
and  tho  change  in  volume  per  unit  original  volume.  Therefore,  the  work  don o  in  compressing 
the  bubble  is  more  than  tho  work  done  by  tho  bubble  in  expanding;  this  difference  in  the  work 
done  represents  a  not  flow  of  heat  energy  into  the  liquid.  When  an  error  in  the  derivation  by 
Pfriem  for  this  thermal  damping  is  corrected,  the  curve  for  the  theoretical  thermal  damping 
constant  at  rosonance  agrees  oxactly  with  Willis’  theoretical  curve  as  given  in  tho  report  by 
Spitzor. 

Pulsating  bubbles  expend  a  portion  of  their  energy  in  tho  form  of  spherical  sound  waves. 
Tho  radiation  damping  is  just  this  loss  of  energy  in  tho  form  of  sound  waves. 

Tho  effect  of  viscosity  on  pulsating  bubbles  in  an  incompressible,  viscous  liquid  is 
understood  through  tho  stress  equations  and  the  boundary  conditions,  rathor  than  tho  Navior- 
Stokos  equation  of  motion.  At  the  bubble  surface,  there  are  viscous  forces  acting  which  exert 
an  excess  prossuro;  this  results  in  the  dissipation  of  onergy. 

Experimental  results  verify  that  the  damping  at  resonance  is  due  to  thermal  arid  radia¬ 
tion,  and  possibly  viscous  damping.  The  discrepancies  botwoen  theory  and  experimental  re¬ 
sults  found  in  measurements  by  Moyer  and  Tamm,  and  Carstensen  and  Foldy  are  duo  to  the 
particular  conditions  of  ti  o  oxporimonts.  The  small  discrepancy  botwoon  the  theory  and  the 
results  of  Bauer  nnd  Lnuor  may  be  due  to  an  incroaso  in  viscous  damping  caused  by  the  non- 
spherical  shapo  of  the  gas  bubble,  and  tho  possibility  that  there  is  some  coupling  between 
the  radial  motion  and  the  shape  oscillations  which  may  remove  some  of  tho  onorgy  associated 
with  the  puisation. 
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